Impurity effects on the spin excitation spectra in a d-wave superconductor 
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The effects of nonmagnetic impurity on the spin excitation spectra in a d x 2_ y 2-w&ve supercon- 
ductor are examined, using the self-consistent t-matrix approximation. It is shown that the impurity 
self-energy acts to shift the position of the resonance peak to low frequencies and broaden the peak. 
While the impurity vertex correction causes a broad spectral weight in the spin gap at the impurity 
concentrations where no clear resonance peak is observed. The gaplike feature still remains in low 
frequency region upon the introduction of impurities. Incorporating these two effects, we find that 
the result is in qualitative agreement with experiments on YaBa2(Cui_ :E Zn a ;)306+i/. 
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I. INTRODUCTION 



The spin excitation spectra Imx of high-T c superconductors is extensively studied by inelastic neutron scatter- 
ing(INS) and a consistent picture has emerged in YBa 2 Cu306+a;(YBCO). A remarkable feature in both underdoped 
and highly doped YBCO is that a sharp neutron resonance peak was observed in the superconducting state(SC) at a 
2D wave vector Q=(7r,7r)[l-4]. Also, in the SC state, Imx is restricted to a small energy range limited by a doping 
dependent energy gapfspin gap) in low frequencies and has sinusoidal dependence on q z , the wave vector perpendicular 
to the Cu02 planes Hjq]. Both the resonance peak and the spin gap disappear in the normal state, and the resonance 
energy E r is found to increases monotonically with the superconducting transition temperature T c |3j,|J, therefore 
they appear to be correlated to the superconductivity. 

A number of theories has been proposed to account for this magnetic resonance. Beyond different modifications, one 
may basically divide the explanations into two classes. First, it may result from the spin-flip quasiparticle scatterings 
across the SC gap, causing an enhancement of the electronic spin susceptibility at a specific energy which compensates 
for the loss of spectral weight below the gap. Second, it may be a consequence of a collective mode in the particle- 
particle channel which couples to neutron scattering through the particle-hole (p — h) mixing in the d-wave SC state ||. 
More particularly, the first class includes i) a BCS gap function wth strong Coulomb correlations |^-^| and a non-BCS 
gap function resulting from the interlayer pair tunneling theory of high-T c superconductivity Q , in the framework of 
a d-wave pairing model, ii) a s-wave order parameter with opposite signs in bounding and antibounding bands formed 
within the Cu02 bilayer [Q. 

Experimentally, it was shown that the superconducting properties are modified by nonmagnetic impurities, espe- 
cially T c is rapidly suppressed by the substitution of the copper ions by zinc ions JD| . A possible interpretation of 
these results is made in term of a d x 2_ y 2 order parameter affected by nonmagnetic scattering, assuming that Zn acts 
as a strong resonant scatterer In this case, nonmagnetic impurities have a strong pair breaking effect p4|JT5[ 

and will modify the spin excitation spectra observed in SC state. 

The purpose of this paper is to study the modifications of the resonance peak and the spin gap upon the introduction 
of nonmagnetic impurites in a BCS d x 2_ y 2 superconductor. We treat the impurities in the dilute limit using the self- 
consistent i-matrix approximation |l6| . Both the impurity self-energy effects and the impurity vertex corrections are 
considered in our calculations. In the pure system, a sharp magnetic peak is observed and the spectrum is limited 
by the spin gap in low frequencies. When the impurity self-energy corrections are considered, we find that the peak 
is broadened and its position is shifted to lower energies. Meanwhile, the magnitude of the spin gap decreases but 
only a negligible contribution to the spin excitation spectrum is found in the gap due to the impurity self-energy 
corrections. On the other hand, the vertex corrections alone induce a broad spectral weight in the spin gap at the 
impurity concentration where no clear reasonance peak is observed, and have only slight effect on both the peak and 
the magnitude of the spin gap. 

The paper is organized as follows. In Sec. II, we discuss the model and study the self-energy corrections. Sec. Ill 
contains the effect of the vertex corrections. We present the conclusion in Sec. IV 
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II. THE SELF-ENERGY CORRECTION 



To consider the modulation of the spin susceptibility along the c axis, we investigate a bilayer system with a interlayer 
hopping t± and the same d x 2_ y 2 order parameter in two layers. The effects of antiferromagnetic correlations within 
and between the layers are considered in the random-phase approximation(RPA) form of the susceptibility. The 
Nambu Green's function of single particles for the pure system in the SC state is given by, 

96 (i) , (1) 

where Ui{do — 1) is the Pauli matrices, i — a or b expresses the antibonding or bonding band. For the quasiparticle 
dispersion, we use £i a ^ = — 2i(cosfc c + cosk y ) — At cosk x cosk y — 2t [cos(2A; x ) + cos(2fcj,)] ± tx — with t ft — 
— 0.2, t jt = 0.25, tx/t = 0.44, /x/t = —1.11 corresponding to a fit to the angle-resolved photoemission data on the 
optimal doping YBCO as used before 0. The order parameter is chosen as Ak = Ao(cosfea; — cosfcj,)/2, where 
Aq = 4T c o and T c o the SC transition temperature. 

The nonmagnetic impurities are modeled by a zero-range potential V and its scattering is treated in the self- 
consistent i-matrix approximation |l6f| . In this approach, two parameters are introduced to describe the scattering 
process: c = l/(7riVoV) and T = Tii/irNo, where Nq and rii are respectively the density of states at the Fermi level 
and the impurity concentration. The impurity-average Nambu Green's function g(k, iu> m ) for single particles can be 
written formally as |l4| , |l6| ], 

#(k,^„) = ^° + Afc " 1+ f!" 3 - (2) 
(iuin ) 2 - A| - ) 2 

The tilde symbol represents inclusion of the impurity self-energy corrections, 

4 l) - u n |« = ^ + (3) 

where we have used the fact that the off-diagonal self-energy vanishes for a d x 2_ y 2 symmetry of the gap function. 
In the single-site approximation, the self-energy is given by = FT- . The impurity-scattering t-matrix T- can 



be calculated from fi""i].[~i"i| . 



T (i) = G { q\uj) r(0 = c 



with G$\w) = (l/7riVo)X;/c Tr [o'o5 (i) ( k ^)]- Th e following calculations are carried out in the unitary limit, c = 0, so 
only the Eo contribution remains. The order parameter A(r, 0) and the SC transition temperature T c in the presence 
of impurities are determined from the gap equation. In the weak-coupling limit, it has been shown that A(r, 0) /Ao and 
T c /T c o draw almost the same curve as a function of T, i.e., A(r,0)/A rj T c /T c0 |13T| . The temperature-dependence 
of A(r, T) is taken to be, 



A(r, T) = A(r, 0) tanh(2V(T c /T) - 1), (5) 
where T c is given by the Abrikosov-Gor'kov formula |jl7j| , 

with ip( x ) the digamma function. 

The spin susceptibility for Matsubara frequencies is calculated from, 



(q,iw m ) = -T2^Tt[i5W(k,iw B )^)(k + q J iw TO + iw„)]. (7) 



Its analytic continuation to the real frequency, giving %q (q, w), is performed using Pade approximants ]l8| . When 
g^> is replaced by , Eq.(7) gives the result for the pure system. The antiferromagnetic correlations in the plane 
J|| and between the planes Jx would renormalize Xo^ • This effect is considered in RPA approximation Jl9| , 
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X W>(q,o,)= . (8) 

with J + (q) = J(q) — J± and J(q) = J||(cosg x + cosq y ). We note that the susceptibility described above comes 
from the p — h excitations of quasiparticles within and between the bonding and antibonding bands. However, the 
susceptibility x p?l (q ! _a£) observed in the neutron scattering is related to the excitations of quasiparticles within and 
between the layers |20| . The relation between them can be obtained using the transformation matrix between the 
states in the layer and band representations. It gives, 

X (U) ^ X (22) - \[x + + X~], X (12) = X (21) * = ^[X + - X-}, (9) 

where d is the distance between two layers, x + — (q> w ) + X (Qj^O an d X~ = X (q> w ) + X (q> lu )- Then 
we have, 

X **(q, u>) = X ™ + X (12) + X (21) + X (22) = X + cos 2 ^+x~ sin 2 ^. (10) 

Eq.(10) implies that the experimentally observed sm 2 (q z d/2) modulation of the INS comes from the transitions of 
quasiparticles between the respective bands. 

In following evaluations, the summation over k and n are performed by dividing the Brillouin zone into 1024x1024 
lattices and by summing from n = —100 to n = 100 in Matsubara frequency uj n = ttT(2ti — 1), respectively. The 
number of input points in Pade approximant is chosen to be 100 and J + (Q) to be 0.85 in unit of t(We will use this 
unit in the following). In addition, we take T c q = 0.1 and T = 0.1T c o. 

Results for Imx p ' l (Q, w) versus to are shown in Fig.l. The continuous line corresponds to the pure system which 
reproduces the observed INS features in the SC state. The dashed (dotted) lines are results with the impurity 
self-energy corrections. To understand the impurity effect, let us first address the origin of the peak for the pure 
system within the d-wave BCS framework, which has been studied in Ref. 0. For a qualitative statement, let 
T = 0, and set the coherence factor to unity, then one has Imxo (Q,w) = ~K^2, k S{ui — Ejp — E^q). The energy 

£to)(k) = E { k l) + E^lg which is the function of the 2D wave vector k has a minimum at E^(k) « 2A = 0.8, 
corresponding to both k and k + Q near the crossings of the Fermi surface and the magnetic Brillouin zone. At the 
minimum Inr^ 1 "''' (Q, u>) has a step and correspondingly Rex^(Q,a») has a logarithmic singularity. In the realistic 
calculations, this divergence exhibits a maximum as shown in Fig. 2 and causes a resonant peak due to the RPA 
renormalization Eq.(8). Meanwhile, there is a saddle point at (0, 7r) in the quasiparticle dispersion, and it leads 
to a logarithmic divergence in ImXo*"^(Q,a;). It arises from the transitions between the occupied states located at 

(0, 7r) and empty states above the SC gap, thus the peak position locates at Esp — A + y A 2 , + (£^jy) 2 - F° r the 
dispersion of quasiparticles considered here, the van Hove singularity of the antibonding band at (0, n) lies at an 
energy = —0.25 relative to the Fermi level and that of the bonding band is £^ = —1.12 due to splitting of the 
two bands. It gives E^ « 0.87, which is close to the energy where Rexo^ (q> w ) is divergent, therefore enhances the 
peak. In fact, these two effects are indistinguishable in the calculations and exhibits only one peak in Itiixq^ (Q, ui) as 
can be seen in Fig. 2. Now, we turn to the impurity effects on the spin excitation spectra. In a ci-wave superconductor 
with resonant nonmagnetic impurity scattering, the SC gap Ao is suppressed [ [l2| , [l3[ and causes the shift of the peak, 
which is basically equal to 2A as discussed above, to low frequencies. Meanwhile, the impurity scattering causes the 
decays of quasiparticle states and leads to the damping of spin excitaitons associated with Ithxq^ (Q, ui). It gives rise 
to the broadening of the peak. Exactly this bahavior is observed in Fig.l. Also one can see from Fig. 2 that the peak in 

haxo^ (Q, disappears gradually upon the introduction of impurities. It is because the impurity scattering will wash 
out the van Hove singularity. Consequently, no clear reasonance peak is observed at large impurity concentrations 
(e.g.r/Ao = 0.08) due to this effect and the damping of spin excitations. Another feature in Fig.l is that no significant 
excitation spectrum weight has been found in the spin gap. The origin of the spin gap in the SC state arises from 
a lack of thermally exciting p — h pairs across the SC gap with transition wavevector Q when the exciting energy is 
lower than the threshold E t h ~ 2Ao. So, though the impurity self-energy produces an increase in the quasiparticle 
scattering rate, it may be not strong enough to cause an observeable enhancement to the p — h excitations across the 
SC gap. We note that the impurity vertex corrections entering ImXo ' (Q> w ) consist of the p — h ladder diagrams 
connected by the impurity scatteing lines, which may allow a strong scatterings and lead to singificant modification 
of the spin gap. 
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III. VERTEX CORRECTION 



In the above calculations, the self-energy from the impurity scattering is considered to include the multiple scattering 
of quasiparticles from the same impurity in the noncrossing manner pl| . Because the dynamical susceptibility 
measured in magnetic neutron scattering is believed here to come from the p — h pair excitations, the multiple 
scattering of particles and holes from the same impurity should be examined. That is, we must include the vertex 
corrections due to the impurity scattering, which is displayed diagrammatically in Fig. 3. The single and double 
arrowed solid lines in Fig. 3 stand for the normal and pairing Green's functions of particles and holes renormalized by 
the impurity self-energy. The dashed line is the impurity interaction and the impurity is represented by a cross x . 
The multiple scatterings in the form of ladder diagrams and the multiple scatterings of quasiparticles from the same 
impurity can be explicitly seen from Fig. 3 (b) and (c), respectively. The vertex-corrected spin susceptibility can be 
written as a 4 x 4 matrix equation plj , 

■ \ M (l: >\q,iuj m ,iu} n ) 

Y„ (q.lUJm = i > , (11) 

° Y l-/W)(^m,^„)MW)(q,zu; m) ia; n )' 

where T(iu> m , iuj n ) = 1 — I^(iu> m , iuj n iui m , iu> n ) is the dressed vertex and the impurity-scattering lines are 

given by, 

I^\lLU m ,iLO n ) = --J£—T£\iu m + lLU n )T^\lLO n ). (12) 

The spin-triplet particle-particle channel into the p — h bubbles by transfroming e.g. a spin-down particle into a 
spin-up hole and vice versa via the mixing with the SC condensate is not included here, because its contributions 
to RPA normalized spin susceptibility is zero when one considers the AF correlations in the form of that in t — J 
model |22j. Thus, the components of M are, 

K«')r„ ,v, i - _ ,v, _i_o„v, _*v, i __j\>(y)/ 



(2tt) 



M$ (q, -iu m , -tw n ) = - / AG^ (p + q, iuj m + i0J n )G^ (q, iw n ), (13) 



M^P (q, iuj m , iuj n ) = M^ j) (q, iw m , iu)„) = - J ^2^ F(l) (p + ^ iuj m + iu n )F^ (q, iuj n ), (14) 

where G' ! '(q, iw n ) and F' l '(<i, iu> n ) ar e the normal and paring Green's functions of quasiparticles which has been 
renormalized by the impurity self-energy. 

The equations (11), (13) and (14) are calculated by using the same method described in Sec. II. Results for 
Imx ph (Q,oj) are shown in Fig. 4 for the same impurity concentrations as those in Fig.l. In contrast to the effect 
of self-energy corrections, an apparent contribution to spin excitation spectra is observed in the spin gap at large 
impurity concentrations where no clear reasonance peak is observed. In order to separate the contribution of vertex 
corrections from self-energy corrections, we have carried out the similar calculations in which the Green's functions 
of the impurity-free system in are used. The result shows that the signal in the spin gap is solely due to 

the vertex corrections, meanwhile the magnitude of the spin gap as well as the position and the width of the peak 
remain unchanged, except for a slight enhancement of the peak height in the presence of only vertex corrections. 
The broad contribution in the spin gap may be understood as the strong scattering involved in the impurity vertex 
which allows for the multiple scatterings due to ladder diagrams. To address the reason why this strong scattering 
takes effect mainly in low frequencies, we show in Fig. 5 the decay rates of the quasiparticles due to the impurity 
self-energy implied by l/r|^ (w) = — 2ImS^(o;) = — 2TT^ l \uj). The similar result has been obtained by Quinland 
and Scalapino |2l]] . We can see that the decay rates increase as the frequency decreases and reache its maximum at 

lu = 0. Because the impurity-scattering lines in the vertex corrections is directly related to Eg 4 ' (oS) as expressed in 
Eq.(12), this enhancement is amplified by the multiple scatterings in the form of ladder diagrams. We note that, in 
the absence of the vertex corrections, this enhancement is not strong enough to lead to an apparent spectral weight 
in the gap as discussed in Sec. II. From Fig. 4, we can also see that the gaplike region in low frequencies still retains 
in the impurity-doping system . We may ascribe it to the fact that the off-diagonal impurity self-energy Si vanishes 
identically for a d x 2_ y 2 order parameter and therefore the angular (e.g. nodal) structure of the SC gap is not changed. 
According to these features, we find that the overall modifications of the spin excitation spectra upon the doping of 
nonmagnetic impurity are in qualitatively consistent with the INS measurement on YBa2(Cui_ y Zn a )3C , 6+; / p3| , p4| . 
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However, the spectral weight in the spin gap is still not large enough to account for quantitatively the experimental 
result |23| . We note that a nonmagnetic impurity such as Zn in the CuC>2 planes is believed to induce a local magnetic 
moment and lead to additional spin-flip scattering p3j . From the above discussion, we think that this scattering may 
lead to more significant spectral weight in the spin gap than that given here. Nevertheless, a detailed investigation of 
this effect is required and will be carried out in future. 



IV. CONCLUSION 



We have calculated the spin excitation spectra below T c for a model d x i_ y i -wave superconductor with resonant 
impurity scattering. The impurity self-energy corrections shift the position of the resonance peak to low frequencies 
and broaden the peak. As the impurity concentration increases, the resonance peak disappears gradually. When no 
clear reasonance peak is observed, the impurity vertex corrections cause a broad contribution to the excitation spectra 
in the spin gap, but the memory of the spin gap still retains. Thus, impurity-scatterings, together with the vertex 
corrections, account for qualitatively the experimental measurement on Zn-doping YBa2Cu306+a;. 



ACKNOWLEDGMENTS 



One of the authors (J.X.Li) acknowledge the support by National Nature Science Foundation of China. 



[1] J. Rossat-Mignod et al, Physica C 185-189, 86 (1991). 

[2] H.A. Mook et al Phys. Rev. Lett., 70, 3490 (1993); P. Dai et al, Phys. Rev. Lett. 77, 5425 (1996). 

[3] H. F. Fong et al, Phys. Rev. Lett. 75, 316 (1995); H.F.Fong et al, ibid. 78, 713 (1997). 

[4] L.P. Regnault et al, Europhys. Lett. 38, 313(1997). 

[5] E. Dernier, and S.C. Zhang, Phys. Rev. Lett., 75, 4126 (1995). 

[6] K.Maki and H.Won, . 72, 1758 (1994); PMonthoux and D.J.Scalapino, Phys. Rev. Lett. 72, 1874 (1994); F. Onufrieva 
and J.Rossat-Mignod, Phys. Rev. B, 52, 7572 (1995); H.Fukuyama, H.Kohno and T.Tanamoto, J. Low Temp. Phys. 95, 
309 (1994). 

[7] G. Blumberg, B.P.Stojkovic and M.V.Klein, Phys. Rev. B, 52, 15741 (1995); D.Z. Liu, Y. Zha and K. Levin, Phys. Rev. 
Lett., 75, 4130 (1995). J. P. Lu, Phys. Rev. Lett, 68, 125 (1992); M. Lavagna and GStemmann, Phys. Rev. B, 49, 4235 
(1994); 

[8] N. Buhit and D.J.Scalapino, Phys. Rev. B, 53, 5149 (1996). 

[9] L.Yin, S.Chakravarty and P.W.Anderson, Phys. Rev. Lett. 78, 3559 (1997). 
[10] LI. Mazin, and V. M. Yakovenko, Phys. Rev. Lett., 75, 4134 (1995). 
[11] G.Xiao et al, Phys. Rev. Lett., 60, 1446 (1988). 
[12] F.Hotta, J. Phys. Soc. Jpn. 62, 274 (1993); 
[13] Y.Sun and K.Maki, Phys. Rev. B51, 6059 (1995). 
[14] P.J.Hirschfeld and N.Goldenfeld, Phys. Rev. B48, 4219 (1993). 
[15] R.Fehrenbacher, Phys. Rev. Lett. 77, 1849 (1996). 

[16] P.J.Hirschfeld, P.Wolfe and D.Einzel, Phys.Rev. B37, 83 (1988); S.Schmitt-Rink, K.Miyake and C.M.Varma, Phys.Rev. 
Lett. 57, 2575 (1986). 

[17] A.A.Abrikosov and L.P.Gor'kov, Zh.Eksp.Teor.Fiz.39, 1781 (1960) [Sov.Phys.JETP 12, 1243 (1961)]. 

[18] H.J.Vidberg and J.W.Serene, J.Low Temp.Phys.29, 179 (1977). 

[19] M.U.Ubbens and P.A.Lee, Phys.Rev.B50, 438 (1994). 

[20] M.Sato et al, Phys.Rev.Lett.61, 1317 (1988). 

[21] S.M.Quinlan and D.J.Scala pino, Phys.Rev.B51 , 497 (1995). 

[22] J.Brinckmann and P.A.Lee, |sond-mat/9710065| ; G.Baskaran and P.W.Anderson, |sond-mat/970607£ . 



[23] Y.Sidis et al, Phys.Rev.B53, 6811 (1996). 
[24] H.Harashina et al, J.Phys.Soc.Jpn.62, 4009 (1993). 

[25] D.Poilblanc, D.J.Scalapino and W.Hanke, Phys.Rev.Lett.72, 884 (1994); A.V.Balatsky, M.I.Salkola and A.Rosengren, 
Phys.Rev.B51, 15547 (1995). 



5 



FIGURE CAPTIONS 



Fig.l. Imaginary parts of the renormalized susceptibility lmx ph (Q,&) versus frequency lu for various impurity 
concentrations r/A . The solid line represents the result of the pure system. Only the self-energy corrections arc 
considered. 

Fig. 2. Imaginary Im^o"^ (Q, oj) and real KeXo"^ (Q, w) parts of the bare susceptibility defined in Eq.(7) versus fre- 
quency u> for various impurity concentrations r/A . a and b represent the antibonding and bonding bands respectively. 
The result for x^ 6) (Q,w) is very similar to Xo"^(Q> w )- 

Fig. 3. Diagrammatic representation of the impurity vertex corrections to the spin susceptibility in the dilute limit 
(see text). 

Fig. 4. Imaginary parts of the renormalized susceptibility Im% p (Q, ui) versus frequency u> for various impurity 
concentrations r/A . The solid line represents the result of the pure system. Both the self-energy and vertex 
corrections are considered. 

Fig. 5. Decay rates of quasiparticles in the antibonding band (w) due to impurity self-energy corrections in 

the unitary limit. Results are shown for various impurity concentrations T/Aq. The result for the bonding band is 
very similar to that shown here. 
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